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^ ! ABSTRACT 

Q 

l/^ ■ A cloud of gas collapsing under gravity will fragment. We present a new theory for 
this process, in which layers shocked gas fragment due to their gravitational instability. 

. Our model explains why angular momentum does not inhibit the collapse process. The 

p ^. theory predicts that the fragmentation process produces objects which are significantly 

Q \ smaller than most stars, implying that accretion onto the fragments plays an essential 

I role in determining the initial masses of stars. This prediction is also consistent with 

^ ' the hypothesis that planets can be produced by gravitational collapse. 



1. Introduction 



Understanding the formation of stars b y gravitational collapse i s a key problern in ast rophysics. 
Sophisticated simulations (see, for example. Bate h Bonnell ( 2005 ): Goodwin et aJ ( 2004 )) are used 



to model star formation within cold molecular clouds in the interstellar medium, but there are sig- 
nificant limitations in the resolution and range of lengthscales achievable with existing computers. 
Furthermore, it is desirable to have simple models which can be used as a framework for under- 
standing the results of simulations as well as fundamental questions, some of which have not yet 
been given satisfactory answers. 

A theory of star formation should explain why stars have a rather narrow distribution of 
masses, despite the very broad distribution of physical properties of the molecular clouds. It should 
also estimate the typical mass of a star and the mass of the smallest objects which are formed 
by gravitational collapse. When a cloud of gas undergoes gravitational collapse, conservation of 
angular momentum implies an increase in its angular velocity. One might expect that this increase 
in angular speed limits the extent of the collapse, and one would like to estimate the importance of 
this effect in limiting the collapse process. Also it is known t hat, on the hig h- mass side of its peak. 



the mass distribution is well approximated by a power-law (|Salpeterjll955l ). an observation which 
theory should explain. 



In this paper we propose a model which gives a simplified but physically well-grounded account 
of gravitational collapse. We show that fragmentation arises because of the gravitational instability 
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of sheets of relatively dense gas produced by shocks. Our theory does not give precise quantitative 
predictions, but it does give a framework which could be used for producing a quantitative under- 
standing of simulations. We find that angular momentum does not create a barrier to gravitational 
collapse. We predict that fragmentation ceases when the fragment size is considerably smaller than 
the mass of a typical star, implying that that mass of a protostar increases by a large factor due 
to accretion of residual gas after the initial formation of a gravitationally collapsed core (and coa- 
lescence of cores could also play a role). There is support from si mulations for the hyp othesis that 
accretion is important for determining the initial mass of stars ([Bate &: Bonnellll2005l ). A theory 
for the initial mass distribution must therefore go beyond the distribution o f mas ses of the gravi- 
tationally condensed cores which we discuss in this paper. iBate &: Bonnelll (|2005l ) have described 
a simple model for this accretion process which gives results consistent with their simulations, but 
a physically well-motivated model will require substantial additional work. 



2. Other approaches to gravitational fragmentation 



Before describing our own approach we discuss some of the existing approaches to gravitational 
collapse. A fundamentally important concept i n most discuss ions of gravitational collapse is the 
Jeans mass (which is discussed and extended by iBonnod (119571 )). A region of initial density po will 
collapse by gravitational self attraction in a time of order the Jeans time 



tcol ~ 



(1) 



provided gravity can overcome the kinetic energy of relative motion. In an initially near-homogeneous 
region, Jeans argued that the size of the gravitational collapsing region is determined by the dis- 
tance that a sound wave, with velocity Cg, can travel in time tcol- This leads to the Jeans estimate 
for the mass produced by gravitational collapse: 

The low gas densities which are ascribed to molecular gas clouds give Jeans masses which are usually 
much larger than the mass of a typical star, indicating that the collapse is accompanied by frag- 
mentation. One theoretical approach to explaining fragmentation stems from the observation that 
if the collapse proceeds isothermally, the J eans rnass d ecreases as the density increases, suggesting 
that the collapsing material may fragment (|Hovldll953l ). It has been argued that the fragmentation 
will cease when the material is no longer able to cool by emitting radiation as it collapses: in the 
case of adiabatic collapse the increase of the speed of sound due to increasing temperature rn eans 
that the Jeans mass ceases to be a decreasing function of density. iLow &: Lvnden-Belll (|l976l ) ar- 
gued that the collapsing cloud starts to behave adiabatically when it becomes opaque at the Wien 
wavelength which corresponds to its temperature. They argued that the minimum Jeans mass is 
insensitive to the opacity of the gas, and used their opacity criterion to estimate the minimum 
mass for a brown dwarf star, finding a result approximately 7 x IO'^Mq, where Mq is the solar 
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mass. The effects of rotation of the collapsing gas are argued to increase this minimum mass, as 
are magnetohydrodynamic effects where these are releyant. 



The arguments put forward bylLow &: Lynden-Belll (1197611 haye been used and ela borated by 
many other authors (for example, Reesl ( 1976 ): Silk ( 1977 ): Masunaga Sz Inutsuka (1999)), but they 
are subject to two major criticisms. Firstly, it can be questioned whether the Jeans mass is a valid 
estimate for the size of the fragments, because the Jeans mass concept is only meaningful when 
the initial density distribution is nearly uniform. This uniformity condition is not satisfied once 
the gravitational collapse has started. A second criticisr n concerns the criterion used to determine 
when the collapsing gas starts to behave adiabatically. iMasunaga &: Inutsukal (|l999l ) argued that 
the criterion of the cloud becoming opaque is not relevant, and that the correct approach is to ask 
when the rate at which thermal energy generated by the gravitational collapse outstrips the ability 
of the cloud to dissipate radiant hea t witho ut an increase in temperature. A further difficulty with 
the approach in lLow &: Lynden-Belll (119761 ) is that it only yields a single estimate for the minimum 
mass of a star, rather than the probability distribution of masses. 

An al ternative picture of gravitetiona l collapse, termed 'turbulent fragmentation' has been 
proposed (jPadoan &: Nordlundl |2002| . 12004 ) . Simulations of gravitational collapse typically show 
the production of many supersonic shock waves, which produce a highly non-uniform density of gas. 
In this approach it is argued that the effects of repeated shocks can produce a log-normal density 
distribution, with regions where the density is orders of magnitude higher than the typical density 
of the molecular cloud. This makes it possible to produce regions where the Jeans mass is orders of 
magnit ude smaller than the solar m ass, and they use their model to predict the prevalence of brown 
dwarfs ([Padoan Nordlundl 12004 ). At large masses, they advance an argument which relates the 
observed power-law distribution of masses to the power-law spectrum of turbulent fluctuations via 
a 'scale-dependent Mach number', but the dynamics of the collapse process is not incorporated 
explicitly. Also, in these models, the critically important parameters of the log-normal distribution 
are derived from ad-hoc assumptions about the properties of the collapsing gas cloud, and they do 
not take detailed account of the role of radiative transfer in controlling the collapse of protostars. 
Our theory also includes the effects of shocks: we give a fuller explanation of the differences between 
our own theories and those of Padoan and Nordlund in the concluding remarks to our pa per (section 
El). Fi nally, here we remark that the turbulent fragmentation model builds upon work by lElmegreen 
(jl997l ). who proposed a theory in which the structure of the collapsing gas cloud is assumed to be 
a fractal set. 

Other discussions of the problem of gravitational collapse emphasise the role of magnetohydro- 
dynamic effects. Magnetic press ure may be very important in stabilising a molecular cloud against 
collapse (jMestel fc: Spitzerill956l ): although the gas is only weakly ionised its effective conductivity 
is high enough to bind the magnetic field lines, and ambipolar d iffusion appears to be too slow to 
allow disconnection of the field lines (IBasu &: Mouschoviaslll994l ). However we argue that once the 
collapse is initiated it is a freefall process in which magnetic and hydrostatic pressures play little 
role until shock waves start to form. 
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3. A new approach to collapse and fragmentation 

Our own theory takes account of the structure of an inhomogeneous collapsing gas cloud by 
specifying a set of points in a multidimensional parameter space. The structure of the fragmented 
cloud can be characterised by considering the set of local maxima of gas density at any instant. 
We assume that gravitational instability will develop each of these density maxima into a discrete 
fragment. Each of these points may be described by the gas density p at that point and also by 
several other parameters. The most important additional parameter is the size scale, L, of the 
density maximum. If a density maximum is associated with a fragment mass M, we define the 
lengthscale by writing M = pL^ . If a density maximum is not well separated from its neighbours, 
we can define its mass by the following procedure. We divide the cloud into Voronoi polygons, 
based upon the set of positions of the density maxima, and determine the total mass M^ew in each 
of these cells. We define the lengthscale L by writing Mf.c\\ = pL^ , where p is the maximum density 
within the cell. 

In addition to p and L, other parameters can be used to characterise the fragments. The most 
important of these is the angular momentum J of the fragment (relative to its own centre of mass), 
which can be used to define a characteristic rotation rate of the cloud fragment through the 
relation J = MLP'Q.. In addition, thermodynamic variables such as the characteristic temperature 
and pressure of the fragment could be used, as well as variables describing its magnetohydrody- 
namic state. However, we are primarily concerned with the behaviour of fragments while they are 
in isothermal collapse, so that temperature is not a relevant variable. Furthermore, at the level of 
approximation which is considered below, there is no significant distinction between a gas which 
dominated by magnetic pressure and one which is dominated by hydrostatic pressure, except that 
in the former case the speed of sound Cs would be replaced by the Alvfen speed. We therefore 
propose to characterise the collapsing gas cloud by a set of points in a three-dimensional space, 
with coordinates p, L, i7, representing the peak density, size scale and rotation rate of the region 
associated with each density maximum. Initially we have just one point, representing the parame- 
ters of the entire molecular cloud. As the collapse proceeds, the cloud fragments and the number of 
representative points increases, until finally each point reaches a region of parameter space where 
fragmentation no longer occurs. This description could be used to characterise empirically the 
results of numerical simulations of gravitational collapse, but here we are concerned with a theo- 
retical description of the dynamics and fragmentation of points in (p, L, $7) space. For simplicity of 
exposition, in the following discussion we initially neglect the role of angular momentum, giving a 
description of the dynamics which is confined to the = plane. We discuss the effects of including 
the rotation rate variable in section [U where we show that the evolution of the rotation rate is 
(approximately) slaved to that of another variable, which characterises the strength of shock waves. 
In this way we show that the effect of centrifugal effects in resisting collapse does not become more 
significant as the collapse proceeds. 

It is convenient to use logarithmic coordinates in (p, L) space, so we shall discuss the motion 
of points in the log p, log L plane, illustrated in figure [TJ As the collapse proceeds, points move 



- 5 - 




Fig. 1. — Fragments of a molecular cloud are described by their density p and lengthscale L. In 
a logarithmic plot of the parameter space, the fragmenting dynamics is confined to a triangular 
region. The three lines bounding this region represent the mass of the molecular cloud, the critical 
density for stability of fragments against evaporation, and the dissipation limit beyond which the 
collapse becomes adiabatic, and fragmentation ceases. A typical trajectory has jumps associated 
with the production of shock waves. 
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downwards (to smaller L) and to the right (to higher p) , and a representative point is replaced by two 
or more points when the density maximum fragments into two or more maxima. If a cloud collapses 
in a homogeneous free-fall, the evolution is along a line of constant mass, M = pL^. Between 
fragmentation events, the representative points therefore move to the right in (log p, log L) space, 
along a straight line of slope — ^, representing evolution at constant mass. When fragmentation 
occurs the number of representative points increases. We shall argue that the dominant mechanism 
for this involves the occurrence of shock waves in the collapsing cloud, which implies a discontinuous 
increase in the density p, accompanied by a discontinuous reduction in the lengthscale L. 

We argue that the dynamical processes leading to fragmentation only occur in a certain trian- 
gular region of the (log p, log L) plane. One of the edges of this triangle is a line of slope — ^ which 
represents a line of constant mass, equal to the original mass of the molecular cloud, Mq = pqLq, 
and all of the points representing fragments of this cloud must lie below this line. Another con- 
straint is that collapse only proceeds if the escape velocity from the fragment, which is of order 
GM/L = y^GpLP, is large compared to the characteristic molecular velocity, which is of the order 
of the speed of sound, Cg. If this condition is not satisfied the fragment can disperse by evaporation. 
This criterion, 

GpL^ 



2 



>K (3) 



(where K is some dimensionless constant of order unity) implies that any fragment which can col- 
lapse gravitationally must lie above a line of slope — ^ in the (log p, log L) plane. The collapse of the 
gas may be resisted by macroscopic turbulent motion with a typical velocity fturbi or by magnetic 
pressure effects, associated with propagation of disturbances at the Alfven speed, va- If either of 
these velocities exceeds the speed of sound, we replace Cg with the largest of these characteristic 
speeds in the following calculations. There is another line, which we term the dissipation line, 
which represents the boundary of the region where the cloud fragments can remain approximately 
isothermal by radia. t ing aw ay heat. Our own criterion for this is different from that proposed in 



Low Sz Lynden-Belll (119761 ): we argue that it is not determined by the condition that the cloud has 
become opaque, and for this reason we use the term dissipation limit rather than opacity limit. 
We will argue that the collapse remains isothermal provided that p^^'^L? remains smaller than 
some characteristic number. We assume that fragmentation ceases once the isothermal condition 
is not met. This implies that the representative points must lie below a line of slope — | in the 
(log p, log L) plane for further fragmentatio n to be possible. We discuss our dissipation criterion 



(which is closely related to one proposed bv lMasunaga k. Inutsukal (|l999l )). and its relation to ear- 



lier work, in section O The fragmentation dynamics is therefore confined to the triangle illustrated 
in figured) once the collapse trajectories exit they no longer fragment. Trajectories which pass the 
dissipation line continue to collapse. If any trajectories pass the evaporation line, the gas cloud 
which they represent cannot collapse further: it may expand slowly or else the gas may be accreted 
by protostars. 



To complete the specification of the dynamics, we must consider the process by which repre- 
sentative points undergo fragmentation, as well as describing the dissipation limit. Both of these 
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points require substantial addition discussion. Before embarking on this, we consider the connec- 
tions between our representation of gravitational collapse and the standard picture based upon the 
Jeans mass. 



In the standard picture of f ragme ntation during gravitational collapse, as proposed by iHoyle 



(|l953l ) and iLow &: Lynden-Belll (|l97d ). attention focusses on just one variable, p, which is used 
to parametrise the state of the collapsing fragments. The lengthscale is always assumed to be 
the Jeans lengthscale, Lj = Cgtcoij where the Jeans collapse time is icoi = (Cp)"^/^, so that the 
dynamics is assumed to follow a line in the plane illustrated in figure [TJ On this line the density 
p and lengthscale L satisfy GpL? ~ Cg, so that the non-evaporation condition is only marginally 
satisfied. In this standard picture, based upon the density dependence of the Jeans mass, the 
dynamics of the collapse process is therefore envisaged to follow the lower edge of our triangle, with 
slope — ^, until the dissipation limit is reached and fragmentation stops. The standard approach 
gives no picture of the mechanism by which fragmentation occurs, nor does it give any framework 
to estimate the probability distribution of fragment sizes, or the role of angular momentum. 

The next section discusses the role of shocks. The dissipation limit for isothermal fragmenta- 
tion is explained in section [5l In section [6] we discuss the role of angular momentum in limiting 
the progress of the fragmentation process. In section [T| we use these results to discuss the prob- 
ability distribution of fragment sizes. A quantitative theory for the fragment sizes would require 
information about the frequency with which shocks occur. The relevant statistics would have to be 
determined by numerical simulation, but it is clear that the distribution of fragment sizes has a 'long 
tail' in the high-mass range, which is similar to the Salpeter distribution. Section [8] summarises 
our conclusions and discusses the comparison between our theory and 'turbulent fragmentation'. 
We argue that the fragment sizes produced by our model are significantly smaller than the sizes of 
typical stars, so that accretion must be invoked to explain the initial mass function of stars. 



4. The role of shocks 



In almost all circumstances the collapsing gas cloud fragments will not be spherically symmet- 
ric, and the first singular structure to appear as the collapse proceeds will be a shock wave. The 
time at which the shock wave appears in collpase of a fragment of typical density p will be of order 
tco\ = ^1 \JGp^ but the exact time at which the shock appears and the its strength will depend 
upon the initial density and velocity distribution of the cloud. The collapse of a spherical cloud of 
uniform density is described by a similiarity solution, in which all of the mass concentrates at the 
centre at the same time, and correspondingly a near-spherical cloud of near uniform density will 
produce a very strong shock at a late stage of the collapse. However the typical case will result in 
a shock strength which is of order unity. In the following we will regard the times and strengths of 
the shock events as random variables with an unknown distribution which we could determine by 
simulation of an ensemble of collapsing gas clouds. 
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Pressure plays little role in the free-fall, until the point at which shock waves are generated. 
There is easily enough energy for freefall to accelerate the gas cloud so that two parts of it will 
meet at supersonic speed, v » Cg. Where this happens the gas at the contact region is compressed, 
and shock waves travel out from the compressed region into each of the incoming masses of gas. 
The shock dissipates most of the relative kinetic energy of the colliding masses of gas into heat. If 
the relative speed of the collision is v, we expect that the thermal molecular velocity in the shocked 
region is comparable to v. Correspondingly, the shocks might be expected to propagate with a 
speed which is some finite fraction of the initial speed v. Note that the collision speed may be 
very large, and correspondingly the temperature of the shocked gas may be very high. 

However, it will be argued that the characteristic radiative cooling timescale, tj-ad; is small 
compared to the Jeans timescale, until the final stages of fragmentation. This indicates that the 
shocked gas loses energy very quickly by radiating heat, so that actually the shock fronts will move 
at speeds which are very small compared to v. In fact, if the gas cools to the ambient temperature 
where the speed of sound is Cg, we expect that the shock front speed Vg is comparable to Cg. The 
collapsing gas therefore forms sheets which have a very high density, relative to the that of the 
initial gas cloud: the density is increased by a factor 

p/p = v/Vs ~ -u/cs . (4) 

The shocked material will build up into a sheet of thickness 

h ~ Cgtcoi (5) 

where tcoi is the timescale over which the gas clouds are in collision. The formation of the shock is 
illustrated in figure [2j 

If magnetic pressure exceeds the hydrostatic pressure, we have to replace Cg by the higher 
speed of magnetosonic waves (the Alfven speed), but the basic principle remains unchanged. That 
is, energy is provided by freefall and efficiently dissipated to an isothermal state by a shock which 
cools rapidly. 

Consider the collapse of a cloud fragment of size Lq and density pQ. As the collapse proceeds, 
freefall will take these values to p, L along a line of constant mass M = pV^. After the fragment 
has collapsed to some fraction / = L/Lq of its original size, a shock wave starts to form. The value 
of / depends on the details of the initial phase-space distribution of the mass within the fragment. 
When the lengthscale is L we expect that the collision velocity of colliding gas clouds is 

V ~ ^GMjL = y^GpL^ . (6) 

The characteristic time over which the shock wave propagates through the gas cloud is the same 
as the characteristic time for gravitational collapse of the fragment, namely tcoi ^ l/\^Gp. Thus, 
the shocks generated by collision of masses of gas when the size of the fragment is L form sheets 
with an area of order and thickness of order 

L' ^h^ tcoiCg ' %= . (7) 
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b 



c 



d 

Fig. 2. — Schematic diagram of a two-dimensional slice through a collapsing fragment. An asym- 
metric fragment (a) collapses in a free-fall until a shock forms (b). The shocked layer subsequently 
puckers due to gravitational instability (c) and fragments (d). The daughter fragments collapse at 
a faster rate than their parent. 
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The shock therefore creates a layer of material of thickness L' and density 

= = (8) 

It will be instructive to express p' / p and L' /L in terms of the ratio M/M^^ where M = pL^ and 
Mj = cl/G^/'^p^/'^ is the Jeans mass corresponding to p. We find: 



p' V _ GpL^ _ G^^p^/^pL^y/^ _( M\ 



1/3 



- 



and 



(9) 



(10) 



L ,/GpI? V M 

Having formed sheets of relatively dense gas, with surface density S ~ p' L' , these will now them- 
selves undergo gravitational collapse. The gravitational instability of a sheet of gas with surface 
density S w hich rotates at a uni f orm a ngular velocity J7 is a problem which has been solved (see, 



for example. iBinney Tremaind (|l980l )). and we base our discussion upon its solution. The sheets 
of gas produced by the shock waves will in general be a curved surface, with varying thickness, 
density and in-plane velocity gradient (which may include shear as well as rotation) , but the results 
for a uniformly rotating flat sheet contain the essential physics. The collapse occurs on a timescale 
*coi ~ ■^here the instability rate 7 describing exponential growth of a density perturbation 
with wavenumber k satisfies 

7^ = 27rGS|A:| - vlk^ - 4^^ (11) 

where is a sound speed characterising the propagation of sound waves in the sheet and is 
the local angular velocity of the material in the sheet, about an axis perpendicular to the surface 



(jBinnev &: Tremaindll980l ) . When this expression yields a negative value for 7 , the sheet is stable 



Consider the interpretation of (jlip for the sheet of shocked material in the case where VL = 
0. For disturbances of the sheet which have a wavelength 27r//c which is comparable to or less 
than the sheet thickness, it is necessary to consider the stability problem in three dimensions (so 
that ([TT]) ceases to be applicable). The largest allowed value of k is therefore of order 27r/L', 
where L' is the thickness of the sheet. The inequality ([3|) then implies that the second term is 
negligible for fragments which are far away from the evaporation line in figure [TJ Equation (jlip 
then implies fragmentation: the larger k, that is the smaller the fragments, the more rapid the 
collapse, indicating that the sheet fragments into pieces of size of order L' . This fragmentation has 
a characteristic timescale tfrag ~ 1/ ^/Gp' , which is faster than the timescale for the collapse of the 
parent fragment by a ratio y^v/cs- When is sufficiently large, however, the sheet is stabilised 
against fragmentation: we return to consider this in detail in section O 

Having created smaller fragments of size L' and density p' , these fragments will themselves 
collapse by freefall, and produce additional shocks. Let us consider the how the shock and fragmen- 
tation process is represented on the diagram of the parameter space in figure [TJ We parametrise 
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the fragments by their density and by their distance from the Jeans hne using coordinates 

X = \ogp, y = log(M/Mj) (12) 

respectively. The shock produces jumps in the values of X and Y . From Q, we see that the jump 
in X is equal to AX = ^Y . After the shock, from ([8]) and (llOp we find that the new value of Y is 



fM'\ 

Y = log ^^J= log 



3 



^•.U°-^Y\y -\y ■ (13) 



c; 



When a shock forms the values of the coordinates X and Y therefore jump to 

X' = X + iy, Y' = ]^Y. (14) 

At each shock the number of density maxima (and correspondingly the number of representative 
points) is multiplied by a factor A''. The largest possible value of N is of order {L/L')"^ = exp(2y/3), 
but the value may be much smaller. In particular, if the rotation rate of the sheet is sufficiently 
large, fragmentation may not be possible at all. Between shocks, the representative points move 
along a line of constant slope, representing freefall at constant mass, 

so that the freefall line has slope ^ in the (A, Y) parametrisation. 

To summarise, the dynamics of the collapse process involves a series of jumps in the (log p, log L) 
plane, which are associated with the formation of shocks in collapsing fragments. A typical tra- 
jectory is shown schematically in figure [H Between the shocks, the representative points evolve 
along lines of constant mass, which take them further from the evaporative stability line, where 
M = Mj. The shocks create new fragments closer to this line. The collapse process may be under- 
stood by following the trajectories of representative points in the (X, Y) plane, with coordinates 
defined by equation (jl2p . The motion of the representative points in these coordinates is illustrated 
schematically in figure [3l 

We emphasise that the foregoing is not to be interpreted as a precise description of the very 
complicated dynamics of the gravitational collapse process, but rather a model indicating the order 
of magnitude of the effects of the shock waves. 



5. Condition for gas to remain isothermal 



The large density enhancement caused by the shock waves discussed in section [5] depends upon 
the shocked material radiating heat so rapidly that it remains approximately isothermal. If the 
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shock waves become adiabatic, the density enhancement is or order unity, and shocks become much 
less effective at seeding fragmentation. We assume that fragmentation ceases when the shocks 
become adiabatic. Variou s criteria for the collapsing gas to jjecom e adiabatic have been discussed 



(jLow h Lvnden-Bell 



1976 



Rees 1976 



Masunaga fc Inutsukalll999l ) . Our own approach is very close 



to that of lMasunaga &: Inutsukal (jl999l ). but not close enough that we can simply quote their results. 



We discuss how our arguments compare with those in earlier works at the end of this section. 

First we argue that the timescale i^ad for an element of gas to equilibrate with its surrounding 
by radiative transport is very short compared to the timescale for gravitational collapse, justifying 
the idea that isothermal collapse is possible. A small element of gas, with volume Ay which is at 
temperature T and mass density p contains internal energy 

'ipkT 



AU 



2mo 



-AV 



(16) 



where mo is the molecular mass, and k is Boltzmann's constant, and we assume that the temperature 
is sufficiently low that the hydrogen molecular rotation spectrum is not excited. This element 
radiates heat to its surroundings at a rate 



(17) 



where a is the Stefan-Boltzmann constant and k is the Planck mean opacity. We see that the time 
constant for the radiative transfer of energy, tj-ad = AU/Q, satisfies 

1 SKuiQaT^ 

trad 3/c 



(18) 



In the molecular clouds of the interstellar medium, the optical properties at the Wien wavelength for 
black-body radiation at the temperature of the molecular cloud (approximately 10 K) are dominated 
by dust particles, and scatt ering is insignificant r elative to absorption. The opacity at T = 10 K is 
estimated (see, for example iBeckwith et a\ ([l990|)) to be approximately 



10" 



2i -1 

m kg 



(19) 



The metallicity of stars varies by about an order of magnitude, indicating that there is at least 
an order of magnitude uncertainty in the value of the opacity. The isothermal speed of sound of 
interstellar gas at 10 K is Cg ~ 190 ms~^ (assuming a mean molecular mass of 2.34). Using these 
values we estimate t^ad ~ 2.3 x 10^*^ s ~ 750 yrs. This timescale is short enough to ensure that 
radiative transfer maintains isothermal conditions until the late stages of the collapse process. 

The rate at which heat is generated in a collapsing cloud is of order 

GM2 



Q 



gen 



Ltcol 



^3/2^5/2^ 



(20) 



The collapse ceases to be isothermal when this rate of heat production exceeds the rate at which 
the cloud fragment can radiate heat. This latter question depends upon whether the cloud is 
transparent or opaque at the Wien wavelength corresponding to the cloud temperature. 
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It turns out that the opticahy thin case is most relevant, and we consider this first. For an 
opticahy thin cloud of volume V = which is at a temperature T + AT and which is surrounded 
by gas at temperature T, the rate at which heat is radiated is 

Qthin = 4.Kpa[{T + Ar)4 - r4]L3 (21) 

where k is the opacity, provided that the optical thickness is small, KpL ^ 1. The collapse ceases 
to be isothermal when AT/T is of order unity. Combining ()20p and (j2ip . the condition for near- 
isothermal collapse to cease is 

provided KpL <^ 1. This condition corresponds to a line of slope —3/4 in the (log L, log p) plane, 
or equivalently to a line in the {X,Y) plane (defined by (jl2p ) which also has slope —3/4. 

Now consider the smallest possible fragment mass in the case where the dissipation limit 
is determined by ()22p . which corresponds to the point at which the dissipation line intersects 
the evaporation line, specified by Combining ^ and (^Z2\i we find that the density at the 
intersection point is 

Pmax '■ ^ • (23) 

It is instructive to note that this estimate is the same as that obtained by equating the gravitational 
collapse time, tcoi = (G/?)"^/^ to the radiative transfer time given by (fTSj) . Using the estimate (fT8]l 
in equations ()23p and ([2]), we find the following estimates for the maximum density at which 
fragmentation occurs, pmax, and the corresponding fragment mass, M^^i^: 

/Omax ~ 3.9 X 10-12 kg m'^ M^i„ « 6 X 10^7 kg . (24) 

With these data the assumption that the gas has not become opaque is satisfied: KpL ~ Kp^/'^CgG^^^'^ 
0.05. 

The minimum mass is small compared to that of a star, by about two orders of magnitude. 
Even accounting for uncertainty in the dimensionless prefactors and that most fragments will be 
significantly heavier than the minimum mass, this estimate indicates that the smallest fragments 
produced in gravitational collapse are not protostars, but rather cores upon which material can 
accrete to become protostars. The minimum mass is, however, comparable to that of a gas-giant 
planet, which makes it plausible that gas-giant planets with eccentric orbit s could have been formed 



by gravitational collapse, rather than by the aggregation of dust grains (IRibas Sz Miralda-Escude 



2007 



Wilkinson fc Mehhgll2008l ). 



The opacity may differ from the estimate (jlSp . and studies of stellar atmospheres indicate 
that the metallicity of stars varies over approximately an order of magnitude. If the opacity k is 
multiplied by a factor 9, then pmax is multiplied by 9^ and the minimum mass of a gravitationally 
condensed core is multiplied by 9~^: that is, higher amounts of dust will favour smaller cores. 
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When the gas becomes optically thick, the transport of heat by radiation is then determined 
by diffusion, and heat may be able to escape from the fragment at a sufficient rate to allow it to 
remain isothermal. We found above that for the value of opacity quoted in (jlSp the dissipation line 
is determined by the optically thin case, but that the extinction length is only marginally larger 
than the cloud fragment. We should therefore also discuss the optically thick case, where diffusive 
transport of radiation is the important mechanism. 

There are two types of radiative diffusion behaviour, depending upon whether the thermal 
radiation is predominantly scattered or absorbed and re-emitted. In the scattering case the radiative 
diffusion constant is Drad ~ ^c, where A is the scattering mean free path, and c is the speed of light. 
However, if absorption dominates the opacity, the diffusion of radiative energy is characterised by 
a diffusion coefficient 



^rad ~ (25) 



1 

a^irad 

where i^d is the timescale for emission of black-body radiation from the absorbing elements, and 
a = up is the inverse extinction length. Using (jlSp . we find the diffusion coefficient 



If the temperature varies slowly on the scale of the optical depth, it satisfies a diffusion equation 
and the radiative heat flux is 

q = ^D,^VT . (27) 
mo 

The the rate of heat production Qgen in a collapsing fragment of size L is of order GM^/Licoh 
where M = pL^ is the mass of the fragment and tcol = (Gp)"^/^ is its characteristic collapse time. 
The interior of the fragment will exceed the ambient temperature by an amount AT. Estimating 
the heat flux as \q\ ~ Q/L?', we estimate AT ~ Quiq/ pkLDj-^d- The condition for approximate 
isothermality, AT/T <^ 1, can then be expressed in the form 

AT Q mo L"^ GpL^ 



T L pkTD^^^ T'radicol ci 



- « 1 . (28) 



^rad^col 

where v/cs is the shock strength parameter introduced in section HI From (|22p , we find that the 
condition AT/T = 0(1) is satisfied when L equal to some constant. We conclude that if 

the dissipation line is determined by the diffusive condition, then the dissipation line in figure [1] 
has slope — |, instead of — |. 

Finally we comm ent on some earlier w o rks w hich have proposed different criteria for the collapse 



to become adiabatic. iLow Sz Lvnden-Belll (|l976l ) argued that the collapse becomes adiabatic when 
the gas becomes opaque on the scale of the Jeans length. This criterion is most widely referred 
to, but it does not appear to be supported by any calculation. These authors also propose that 
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in the period before the collapse becomes adiabatic, the temperature of the gas rises as T oc p^^^ , 
rather than remaining isothermal, wh ich seems t o be due an incorrect application of the energy 
balance relation for a volume element. iReed (|l976l ) gives a discussion of what is essentially our own 
criterion for the dissip ation line, but his paper de sc ribes the results as if they ar e consistent with 
the opacity criterion in lLow &: Lvnden-Belll (jl976l ). iMasunaga &: Inutsukal (|l999l l were the first to 
point out that the opacity condition for adiabaticity is not correct, and they appear to be the first 
to point out that diffusion may limit the rate of radiation transfer when the gas cloud is opaque. 
Their discussion assumes implicitly that the collapse process is parametrised by a single variable 
(the gas density). Our own results extend theirs by describing the dissipation limit as a line in the 
(log p, log L) plane. 



6. The role of angular momentum 

If a cloud of gas is rotating, the centrifugal effect resists gravitational collapse. The cloud 
can only collapse if the gravitational collapse rate \/Gp is sufficiently large relative to the angular 
velocity of the cloud, Thus collapse only proceeds if 

% > A" (29) 
where K' is a dimensionless constant of order unity. 

The angular momentum J of an isolated gas cloud is conserved, suggesting that as the cloud 
collapses the rotation rate ^1 will increase and eventually the inequality (j29p will not be satisfied, 
so that gravitational collapse will cease. If the size scale of the cloud is L and the density is p, then 
the angular momentum J is of order J ~ ML'^Q, where M ~ pL^ is the mass of the cloud. In a 
freefall collapse of the cloud, the conservation of angular momentum implies that oc L~'^, and the 
conservation of mass implies L oc p~^^^, so that O oc p^/^. It follows that Gp/i}'^ oc p~^^^, implying 
that the freefall collapse of a gas cloud will be halted when the density p becomes sufficiently large. 

If the gas cloud undergoes fragmentation, however, centrifugal effects can become less of a 
barrier to the collapse proceeding. Note that when we characterise a fragment of a gas cloud, we 
should be concerned with its angular momentum about its own centre of mass. Thus if a cloud 
with angular momentum J breaks up into fragments with angular momenta Jj, the principle of 
conservation of angular momentum does not imply that ^ ■ J, = J. This is because each of the Jj 
is measured relative to an origin which is different from that used to determine J, so that in fact 
we have J > Jj. 

When a cloud fragments by gravitational self-attraction, we expect that as the fragments start 
to become distinct from the original cloud they will have approximately the same angular velocity 
as that part of the cloud from which they form. Although there is no exact conservation principle 
which is applicable, we can assume that the rotation rate of each fragment, Oj, is comparable to 
the rotation rate of the parent cloud: ~ il.. 
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The fragmentation mechanism described in section HI which involves the production of strong 
shocks, produces fragments which have a greater density than the parent cloud, but the rotation 
rate of the gas will not be significantly changed by the formation of a shock. Equation ()29p implies 
that the increased density favours gravitational collapse. 

There are thus two competing effects. A freefall collapse inhibits gravitational collapse by 
increasing the centrifugal effect, whereas the shock-induced fragmentation mechanism favours grav- 
itational collapse by increasing the density. We now consider how to quantify the balance between 
these opposing effects. To this end we introduce a further logarithmic variable characterising the 
state of a gas cloud or one of its fragments: 

z = (f ) ■ m 

The cloud can collapse if Z exceeds a threshold value, which is of order unity. We argued that in 
freefall, Gp /il^ oc p~^^^, so that in the freefall phases the relation between Z and X (defined by 
l^)is Z = Zo-^X (where Zq is a constant determined by the initial conditions). In freefall, we 
also have the relation Y = Yq + ^X (see equation (jlSp : Yq is another constant), so that if a cloud 
undergoes freefall such that X changes by AX, the corresponding changes in Y and Z are related 
by 

AZ = --AY = --AX . (31) 
3 3 

When the shock and fragmentation process described in section S] occurs, however, the density 
increases with remaining constant, resulting in an equal jump of both X and Z, specified by 
equation (jl4p . so that the value of Z after the shock is 

Z' = Z + -Y . (32) 
3 

Equations (j3ip and (j32p show that the evolution of the variable Z, which describes the effect of 
centrifugal forces, is slaved to the evolution of Y, which parametrises the strength of shocks. In 
the shock, equation (jl4p also shows that the value of Y is halved. Thus we see that the evolution 
of Z is coupled to that of Y. During periods of freefall, X changes by increments AXi, and there 
is a corresponding change AYi = ^AXi in the value of Y. Let the values of Y just before shock 
events be Yi. The change in value of Y over many shocks and freefalls is 

y(Xf ) - Yixo ~ Yl l^^^ - ■ (33) 

i 

The value of Y is expected to reach a stationary distribution, in which case ^(-'^f ) — ^(^i) remains 
bounded, so that ([33]) implies that 

{AXi) = {Y,) . (34) 

By analogy with (j33p . the change in Z is 

Z{X() - Z{XO ~ -l^^^ + l^i (35) 
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Using (|34p. we conclude that the difference between the final and initial values of Z has no drift. 
This is consistent with Z also approaching a stationary distribution. 

Although centrifugal effects may stop some fragments from undergoing further fragmentation, 
our results are consistent with the hypothesis that angular momentum does not become more 
significant as the collapse process progresses. 

7. Statistical model for masses 

An advantage of considering the dynamics of the collapse process in greater detail is that it gives 
a framework within which the distribution of fragment masses can be predicted. A quantitative 
calculation of the fragment mass distribution from our model is a difficult problem which will require 
substantial additional effort. However it is of interest to consider how this calculation can proceed 
in principle. In the following we discuss in outline how the fragment mass distribution could be 
calculated. Our discussion indicates that in the limit where the size of the initial molecular cloud 
becomes very large, the distribution of fragment masses becomes universal. We also argue that our 
approach is consistent with a power-law for the large-mass tail of the distribution, analogous to the 
Salpeter distribution (although the typical fragment mass is much smaller than the typical initial 
stellar mass). 

The probability distribution of masses of the fragmentation products is determined by the 
density of points in the (X, y) plane at the point where the fragmentation ceases, that is when 
the representative point passes through the dissipation line. Given this density, w{X^Y\ we can 
determine the probability distribution Pf(lf) for the final value of y, that is the value If where 
the fragment trajectory passes the dissipation line. This distribution can then be transformed to 
determine the final fragment masses, M = Mjexp(y). We first consider the form of the density 
y), and then discuss how this can be used to determine the fragment mass distribution. 

We can define a stochastic process to model the variable Y{X\ We consider a process where 
y(X) grows linearly with X (with gradient e, say) until a shock occurs, at coordinate Xj. Let the 
probability of an interval of growth between shocks being greater than AX be P(AX). When a 
shock occurs, the value of Y is multiplied by a fraction ^. We have seen that e = ^ and [i = \ 
are the appropriate values to use with our model. If the molecular cloud is orders of magnitude 
more massive than a single star, there will be many fragmentation events during the gravitational 
collapse, and in this case the steady-state probability distribution of our stochastic process will 
approach a limit, with probability density function p(y). The number of fragments is expected 
to increase exponentially as the fragmentation proceeds, so that the density of fragments in the 
(X, y) plane which results from fragmentation of a very large body may be approximated by 

w(X, Y) = exp(r/X)p(y) (36) 

where ry is a universal constant and p is a universal function. 
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Now consider how the fragment mass distribution may be obtained from (j36p . We argued that 
the dissipation hmit is determined by the relation = const., which corresponds to a hne of 

slope — I in the {X, Y) plane, with coordinates defined by equation (|12p . Fragmentation stops when 
the representative point describing a fragment crosses the dissipation line, at the point (Xf,Yf). 
We want to compare the mass Mp of this fragment with that of the minimum mass fragment, Mmin, 
corresponding to the rightmost apex of the triangle in figure [H at coordinates (^max)O). Because 
the slope of the dissipation line in the (X, Y) coordinates in — |, we have X{ = Xmax — |^f- Noting 
that the Jeans mass Mj is proportional to p~^^'^ , we have 

exp(Yf) = exp exp(yf) 



= exp (^-Y^ . (37) 

To determine the distribution of fragmentation product masses we must therefore evaluate the 
distribution of the coordinate Yf and use equation (j37p to relate If to the fragment mass. The 
density at the position (Xf,lf) relative to that at the corner of the dynamical triangle is: 

- exp --r/Yf — — . (38) 



lt^(Xmax,0) ^ V 3 ' 7 p(0) 

It is a difficult task to determine the probability P(AX). In the following we consider the 
result of applying the simplest model, which is to assume that the probability of going from X to 
X + 5X without the formation of a shock is independent of X (that is, we assume that AX has 
a Poisson distribution). In this case P(AX) = exp(— AAX), for some constant A. We are really 
interested in the distribution for the case where fx = ^, but it is difficult to achieve a solution in 
that case. Rather, we discuss the behaviour of the solution in the limit as ^ 0, which is very 
easily obtained (as we show in the following paragraph). Numerical studies show that when P(AX) 
is Poisson distributed, the tail of the distribution is very similar when /j = ^. 

In the limit where // = 0, the coordinate Y is reduced to zero every time a shock occurs. In 
this case the probability that Y is less than Yq is 1 — P(Yq/€). Let us consider the distribution of 
fragment masses in the case where P{AX) is approximated by a Poisson distribution, so that the 
probability density of Y is then p{Y) = Aexp{—XY/e) (where A is a normalisation constant). If Y 
has this exponential distribution, the corresponding distribution for M is a power law. Recalling 
that the fragment mass Mf is related to Y by equation ([37j) , the probability distribution of fragment 
masses is of the form 

PM(Mf)~C — ^ (39) 



in the large-mass limit, where C is a normalisation constant (note that we have set e = ^ here). The 
distribution of separations between random events is often characterised by a Poisson distribution, 
and we have seen that such a distribution gives a power-law tail of the mass distribution, which 
is analogous to the Salpeter distribution for the masses of stars. However, we observed that the 
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masses of the fragments produced by gravitational collapse are much smaller than the masses of 
protostars. The initial accretion of mass onto the fragments generated by gravitational collapse is 
not well understood, and it is uncertain how or whether the distribution of initial masses is reflected 
in the final masses of protostars. Numerical studies of the collapse process are required to identify 
the value of r] and to determine how well p{Y) is modelled by an exponential distribution. 



8. Concluding remarks 



We have described a mechanism whereby a gas cloud fragments due to the formation of shocks. 
When these occur, only a finite fraction of the material in the cloud becomes shocked and undergoes 
fragmentation. The remainder of the gas remains unaffected, and continues to evolve at a slower 
rate than the shocked material. This residual material may itself form further shocks as it collapses, 
or it may become stabilised against further collapse by centrifugal effects. Whatever is the detailed 
sequence of the collapse process, when non-fragmenting dense cores form, they will do so in an 
environment of low density gas, which was not able to take the fastest route to gravitational 
collapse. This reservoir of low density gas can end up feeding the growth of protostars by accretion 
processes. A theoretical description of this late stage is outside the scope of the theory presented 
here: we cannot readily estimate how much the mass of a nascent protostar will be increased by 
accretion processes, nor the effect of the accretion on the form of the initial mass function of stars. 

It is instructive to c ompa. r e our approach with the turbulent fragmentation theory advanced 
by iPadoan Nordlundl (j2002l . |200J). Shock waves play a central role in both theories, but in 
other respects the theories are very different. The turbulent fragmentation theory suggests an 
origin for the Salpeter mass distribution, a power-law distribution for the mass of heavy stars. The 
power-law in the mass distribution is proposed to result from a power-law relation between relative 
velocity and separation ii i the gas clo ud, analogous to that which characterises fully-developed 
turbulence (discussed by iFrischI (119971 )). This implies that the sho ck strength is predicted t o 
increase with the size of the structures in the turbulent gas cloud (jPadoan &: Nordlundl |2004| ) . 
By way of contrast, in our model the shock strengths approach a universal distribution which 
is independent of the size scale. Also, we predict that the initial mass of fragments is much 
smaller than the typical stellar mass, with the protostellar mass determined by subsequent accretion. 
In the turbulent fragmentation theory the peak of the mass distribution is determined by some 
adjustable parameters, whereas in our model the fragment mass is determined by radiation transfer 
considerations. 



Our model is consistent with the numerical simulation bv lBate &: Bonnelll (|2005l ) in that the 
gravitationally bound cores which form protostars are initially much smaller than typical stellar 
masses. One point of difference is that our model predicts that sheets of dense material formed 
by shocks will fragment into clouds of material with a size similar to the thickness of the sheet, 
whereas it is usually reported that simulations show sheets tending to break up into filamentary 
structures. In our discussion of the gravitational instability we assumed that the local motion of 
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the shocked material is purely rotational. We hypothesise that filamentary structures arise when 
the local velocity gradient in the shocked layer is anisotropic. 

It has been suggested that some of the puzzling features of extrasolar planetary systems, 
such as the high probability of finding highly eccentric orbits, can be explained by the hypothesis 
that planets formed by gravitational collap se (IRibas &: Miralda-Escuda 120071 : IWilkinson &: Mehlig 
2008 ) rather than by the standard model ( Safranov 19691 ) involving aggregation of dust particles 
in a circumstellar disc. According to our theory, the size of the smallest fragments produced by 
gravitational collapse are comparable to the sizes of gas-giant planet s, which lends support to the 
picture of planet formation described by IWilkinson &: Mehligj (j2008l ) . In this picture the planets 
would arise from dense cores which were not able to grow by accretion because of competition from 
a larger core nearby, which grows to become a protostar. 
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Fig. 3. — Schematic illustration of the motion of a representative point in the coordinates X = log p, 
Y = log(M/Mj). The blue lines are intervals of free-fall collapse, which terminate with formation 
of a shock. The shocked surface forms fragments which have new coordinates determined by (114p . 
Fragmentation ceases when the trajectory meets the dissipation line at {Xf,Y{), which in these 
coordinates has slope — |. 



